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A PARABOLIC ANALOGUE OF THE HIGHER-ORDER COMPARISON 
THEOREM OF DE SILVA AND SAVIN 

AGNID BANERJEE AND NICOLA GAROFALO 


Abstract. We show that the quotient of two caloric functions which vanish on a portion of 
the lateral boundary of a H k+a domain is H k+a up to the boundary for k > 2. In the case 
k = 1, we show that the quotient is in H 1+a if the domain is assumed to be space-time C 1,a 
regular. This can be thought of as a parabolic analogue of a recent important result in IDS1I . 
and we closely follow the ideas in that paper. We also give counterexamples to the fact that 
analogous results are not true at points on the parabolic boundary which are not on the lateral 
boundary, i.e., points which are at the corner and base of the parabolic boundary. 


1. Introduction 


The classical comparison theorem states that two nonnegative harmonic functions which van¬ 
ish on the boundary of a Lipschitz, or more in general a NTA domain, must vanish at the same 
rate. An important consequence of this result is that the quotient of two such functions is, in 
fact, Holder continuous up to the boundary (only the function in the denominator needs now to 
be nonnegative now). In some recent remarkable works De Silva and Savin have established a 
higher-order version of such result. Specifically, they have proved in [DSl] the following. 


Theorem 1.1. Let D be a C k,a domain in M n , with 0 £ dD. Let u, v be two harmonic functions 
vanishing on dD fl 11(0,1). Furthermore, u > 0 in D and u = 1 at some interior point in D. 
Then, 


( 1 . 1 ) 


C k ’ a ( B { 0,1/2)) 




The classical Schauder estimates imply that u , v are C k,a up to the boundary. Then, by the 
Hopf Lemma we have u u > 0, and from this one can assert that the quotient ^ is C k ~ 1)OL up to 
the boundary. However, Theorem 11.11 remarkably states that the ratio is in fact C k,a up to the 
boundary. The case k = 0 of this result is the boundary Harnack principle mentioned in the 
opening, see jCEMSj and [JK| . 

The purpose of this note is to generalize Theorem 11.11 above to the heat equation and, more 
generally, to linear parabolic equations with variable coefficients. The main results are Theorem 
13.II and Theorem 14. 51 below. Although our work has been strongly motivated by that of De Silva 
and Savin, it has nonetheless required some delicate adaptations to the parabolic setting. 

It is worth mentioning here that, besides being an interesting regularity result in its own 
right, a direct application of Theorem 11.11 above implies C°° smoothness of a priori C l,a free 
boundaries without the use of the hodograph transformation as in [KNj . [K.XSl . a tool that so 
far has been the standard way of establishing smoothness of free boundaries starting from C 1,a . 
For this aspect one should see Corollary 1.2 in [DSl] . Having said this, we would like to mention 
that the hodograph transformation in |KN) . [KNS] does in fact imply real-analyticity of the free 
boundary, which is instead not implied by Theorem ll.il Nevertheless, Theorem 11.11 provides a 
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new perspective in the study of Schauder theory and free boundary problems. Theorem II. II has 
also been extended to slit domains in [DS2j. In the same paper such result has been used to 
establish smoothness of the free boundary in lower-dimensional obstacle problems of Signorini 
type near regular points. The real analyticity of the free boundary near regular points in the 
elliptic thin obstacle problem has been recently established in [ KPS ] by using a method based 
on hodograph transformation. 

These recent results and their applications to free boundary problems motivated us to inves¬ 
tigate their parabolic counterpart. Our main result is Theorem 13.11 below which constitutes the 
heat equation counterpart of Theorem 11.11 We mention that the case k = 0 of Theorem 13.11 can 
be found in [ FSY] . In the present paper we make the observation that the ideas in |DS1| can be 
successfully adapted to the parabolic situation. The idea of the proof in [DSl] is to approximate 
(after a suitable change of coordinates) v by polynomials of the type x n P by a compactness 
argument which uses Schauder estimates. Then, finish by remarkable idea that x n P can in 
fact be replaced by uP. In our situation, as in Theorems 13.11 and 14.51 we show that the ap¬ 
proximating polynomials in the space variable in |DS l] can be suitably replaced by appropriate 
approximating parabolic polynomials, and one can argue in a similar manner. Modulo some 
delicate details, which we have tried to illustrate as much as possible. Similarly to Corollary 
1.2 in [ESI], Theorem 13.11 implies, in particular, the C°° smoothness in the parabolic obstacle 
problem of a C 1,a free boundary near the regular points considered in Theorem 13.1 in [CPS j. 
We note nonetheless that, analogously to the elliptic case, one can establish the space-like real 
analyticity of the free boundary by employing the hodograph transform as in [CPSJ . It remains 
to be seen whether the analogue of Theorem 1.1 in |DS2j is true for parabolic equations since 
such result would have important applications to the parabolic thin obstacle problem which was 
systematically studied in [DGPT] . This question will be addressed in a future study. 

This paper is organized as follows. To better demonstrate the ideas we first establish in 
Section [2] the higher-regularity result in the case k = 1 and for the heat equation. In Section [3l 
still for the heat equation, we analyze the case k> 2. In Section [4] we extend the results of the 
previous sections to non-divergence form operators with variable coefficients. Section [5] closes 
the paper. In it we present an application of Theorems 12.21 and 13.11 to the parabolic obstacle 
problem studied in [ CPS ]. 

Acknowledgment: The paper was finalized during the first author’s stay at the Institut Mittag- 
Leffler during the semester long program Homogenization and Random Phenomenon. The first 
author would like to thank the Institute and the organizers of the program for the kind hospitality 
and the excellent working conditions. 

2. H 1+a REGULARITY FOR THE HEAT EQUATION 

In this section we establish the case k = 1 of the main higher regularity result which is 
Theorem I3T1 below. Our main result is Theorem 12.21 In order to state it we need to introduce 
some preliminary notation and hypothesis. 

With x £ R", t £ R, we will denote by (x, t) a generic point in R n+1 . If xo £ R n and to £ R 
we indicate with 

(2.1) Q r (x 0 ,to) = B r (x 0 ) x (t 0 - r 2 , t 0 ] 

the parabolic cylinder “centered” at (xo,to)- Given an open set G C R n+1 we say that a 
point (xo,to) £ dG belongs to the parabolic boundary of G , and write (xo,to) £ d p G, if for 
every r > 0 the open cylinder B r (x o) x (to — r 2 ,to) contains points of the complement of G 
(notice that (xo,to) 0 B r (x o) x (to — r 2 ,to)). Thus, for instance, when G = Qx (0, T), then 
b p G = (II x {0}) U(dH x [0,T]). We denote by SG the lateral boundary of G, see p. 13 of [Li] 
for the relevant notions. The reader should also see p. 5 of u for the notion of parabolic norms 
and distance. For the for the definitions of C k,a spaces, norm and seminorm, we refer the reader 
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to p. 90 in [GT] . We also refer to p. 46 in [Li] for the relevant notion of H k+a spaces and p. 75 
in [Li] as well for the definition of domains with H k+OL boundaries. 

We now consider a connected bounded open set G C {(x,t) G M n+1 | t < 1}, and we assume 
that (0,0) € d p G. We also suppose that dG D Q 2 (0,0) = SG D Q 2 (0, 0), and that G fl Q 2 (0,0) 
is space-time C 1,a regular, i.e., there exists / G C 1,a such that 

G n q 2 ( 0 ,0) = {(x, t) G < 52 ( 0 , 0 ) I x n > f(x',t)}. 

If we introduce the following notations: 

Ft = {x G R n I (x, t) G G n Q 2 ( 0 ,0)}, Gt = {x G M n | (x, t) G dG n Q 2 (0,0)}, 

then for each t < 0, the set Gt is a (n—l)-dimensional C 1,a submanifold which can be equivalently 
characterized in the following manner 

G t = {(x',x n ) | x n = f(x',t ) and (x,t) G Q 2 (0,0)}. 

For each x G Gt we denote by vt{x) the inward unit normal to Ft at x. In the following discussion, 
whenever there i son ambiguity about the point (x,t) that is being considered, we will simply 
write v instead of vt{x). 

We will use the notation D'f when referring to the gradient of / with respect to the variable 
x' = (xi, ...,x n _i) G M n_1 . Without loss of generality, we may and will assume that 

(2.2) p = p o (0) = e n , /(0,0) = 0, £>7(0,0) = 0, and ||/|| c i,« < c 0 , 

where Co is a dimensional constant which is chosen sufficiently small in such a way that (e n , —3/2) G 
G. This latter property can always be achieved by suitably scaling the domain as we describe 
later. We also normalize the function u appearing in Theorem 12.21 below so that the following 
holds 

(2.3) u(e n , —3/2) = 1. 

Remark 2.1. Hereafter in this paper when we say that a constant is universal we mean that it 
depends only on the dimension n, and the parameters a and k in Theorem \3.1\ below. We notice 
that in the next Theorem 12.21 we are taking k = 1, and thus in this case the dependence would 
be only on n and a. 

The main result of this section is the following H 1+a regularity result. 

Theorem 2.2. Let the domain G C M n+1 be as above and satisfy the assumption (12.21) . Let u 
and v be two solutions to the heat equation in GnQ 2 (0, 0), with u, v vanishing on d p GC\Q 2 (0, 0), 
and suppose that u > 0 in G fl Q 2 (0,0) and that it satisfy the normalization (12.31) . Then, for 
some universal C > 0 one has 

V 

( 2 -4) l|-||// 1 +“(GnQi(o,o)) < C(ll ?; llL° o (GnQ 2 (o,o)) + !)■ 

Remark 2.3. It is worth mentioning here that, although in Theorem, 12.21 we assume that the 
domain be C 1 ’ 01 , instead of just H l+a regular, it is not restrictive in its application to free bound¬ 
ary problems (See Corollary 1,5.31 below). This follows from the fact that the classical boundary 
Harnack inequalities imply that for the parabolic obstacle problem studied in [CPS] the Lipschitz 
free boundary near a regular point as in Theorem 13.1 in fCPSj is shown in the subsequent The¬ 
orem 14.1 in the same paper to be space-time C 1,a regular. The hypothesis of C 1,a regularity in 
Theorem\ 111.HI is only used to apply the Hopf Lemma as in Theorem 3’ in |LNl . Note that the 
H 1+a regularity assumption on the domain does not imply that (4-6) in m holds, which is 
precisely why we assume that the domain be C 1,a regular. 

Remark 2.4. We now illustrate by an example in El that Theorem \2.3\ cannot possibly be 
true at the base or at the corner points of a smooth cylinder. Consider to fix the ideas G = 
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5(0,1) x (0,1). Let u,v be the solutions of the heat equation in G corresponding to Cauchy- 
Dirichlet data g(x,t) = t a and h(x,t ) = tfl respectively. Clearly, u and v vanish at t = 0. 
Assume that f3 < a and denote by K(x,t,y, s) the kernel function for G at the boundary such 
that 


u(x,t)= / K(x, t, y, s)g(y, s)da(y)ds, v(x,t)= / K(x,t,y, s)h(y, s)dcr(y)ds 

Jo JdB Jo Job 

Then, we have 


v(x,t) _ /q Jqb K ( x > Vi s)s f} d<j(y)ds 1 f* f QB K(x, t , y, s)s a dyds _ 1 
u ( x it) /o I dB K (x,t,y,s)s a da(y)ds ~ t a ~h f* f gB K(x,t,y, s)s a dyds t a ~ P 

We conclude that the ratio ^ cannot possibly be bounded as t —» 0 + . This example demonstrates 
that Theorem \2.2\ is not true in a neighborhood of a point (®o>0) € B x {0}. 

The same example can be modified to demonstrate that Theorem \2.2\ is not true in a neigh¬ 
borhood of a corner point (xo,0) € dB x {0}. Let <f> be a smooth function on dB such that <f> 
vanishes in a neighborhood of x o, and let this time g(x,t) = <f>(x)t a , h(x,t ) = (f(x)t 13 . As above, 
we obtain and therefore the ratio ^ is not bounded in a neighborhood of {x o,0). 


Before proving Theorem 12.21 we make some preliminary considerations and reductions, and we 
establish a crucial auxiliary lemma. With u, v as in Theorem 12.21 by Schauder theory (see for 
instance Theorem 4.27 in [O]), we have that u,v € H 1+a up to 8 p G fl Q 3 / 2 (0, 0), say. Moreover, 
by the Hopf lemma in Theorem 3’ in (LN] , the Schauder type estimates, (12.21) . the normalization 
condition (12.31) and the interior Harnack inequality for parabolic equations, we have that 


(2.5) u v (x, t) > c > 0 in d p G D Qi(0, 0), 

where, we recall, v indicates the inward normal at (x,t) € Df. After parabolic dilations of the 
domain, i.e., by considering the rescaled functions 


( 2 . 6 ) 


u ro (x,t ) = 


u(r 0 x,r$t) 


v ro {x,t ) = 


v(r 0 x,r%t) 


r 0 r 0 

we see that u ro , v ro solve the heat equation in the rescaled domain 

G ro = {(x,t) € M n+1 | (tqx, rgt) e G}. 

Moreover, G r ° is given near (0, 0) by the graph of f ro (x',t) = flLhLiLoll. From (12.61) it is easy to 
see that for every (. x , t), (y, s) € G r ° D Q 2 ( 0 ,0) one as 

\D x u ro (x,t) — D x u ro (y, s)\ arn 1 

(| x _„|2 + | i _ s |)./2 £ <■« 

Therefore, given 5 > 0, we can choose the scaling parameter vq = ro([D x u] Ha ^Q n Q 2 ^ 0 0 ^, 5) > 0 
in such a way that 


( 2 -7) [ DxUr 0 \H°‘(G r onQ 2 (0,0 )) < <5- 

It is then clear that for each 6 > 0 the corresponding function u ro satisfying (12.71) does depend 
on 5. We also note that because of (12.21) . if for a given d > 0 the scaling parameter ro is suitably 
chosen, and if u is multiplied by an appropriate constant depending on c in (12.51) . we can ensure 
that the following holds 


(2.8) \\frolic 1 ’* < DxU ro { 0,0) = e n . 

We note in passing that the first inequality in (12.81) represents a flatness assumption of the 
boundary of G which will become important in establishing (12.321) below. Moreover, the choice 
of 6 which is to be fixed later will be determined by Lemma [2.61 where a compactness argument, 
in which we let 5 —> 0, is employed. More precisely, in the proof of Lemma [2.61 5 is chosen small 
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enough so that (12.201) and (12.231) below hold for some p > 0 universal. In conclusion, given the 
choice of 5 determined by Lemma l2.61 from now on to simplify the notation we will let u ro = it, 
v ro = v, G r ° = G , and establish our results for these new u, v and G. 

We now introduce the relevant notion of approximating affine function. 

Definition 2.5. We call P an approximating affine function at (0,0) if it has the form P(x) = 

Ya =i a i x i + °o, with a n = 0. 

As in [DSlI . we have the following intermediate lemma. 

Lemma 2.6. Assume that for some r < 1 and P an approximating affine function with |oj| < 1, 
one has 

( 2 -9) lb ~ “-P||L°°(GnQr(o,o)) <r 2+a . 

Then, there exists an approximating affine function P such that for some C, p > 0 universal, we 
have 

(2-10) H-P - ^||L°°(GnQ 7 .(o,o)) < Cr 1+a , 

and 

(2-11) |b - uP |U°o( G nQ pr (o,o)) ^ ( P r ) 2+a • 

Proof. We let G = {(^, j (. x,t ) € G}, and consider the function v defined in G by the follow¬ 
ing equation 

(2.12) v(x,t) = u(x,t)P(x ) + r 2+a v(—, -^). 

Although in the next definition u does not have the same meaning as v in (I2.12|) above, for later 
purposes we nonetheless abuse the notation and set 

(2.13) ^ 

r 

Before proceeding we pause to recall that in the reduction which precedes Lemma 12.61 given 
any 5 > 0 we have chosen vq > 0, depending on 6, such that the function u = u ro satisfy (12.71) 
and (12.81) . Since we also set v = v ro , it is clear that both u and v do depend on 6, and therefore 
so do u and v. The reader should keep this in mind when below we let 5 —> 0 along a sequence. 
Since it, v and P are solutions of the heat equation, we easily obtain from (12.121) 

n-l 

(2.14) 0 = (Au - v t )(x,t) = 2 ajDjufx, t ) + r a (Av - -j), 

i —1 

where we have denoted by Du = ( D\u ,..., D n u ). Moreover, by (12.81) we know that Diu( 0,0) = 0 
for i = l,..., n — 1, and since by (Hud we also have |Dit]^ a ( G nQ 2 (o,o)) ^ $■> we conclude that for 
all i = 0, ...n — 1, 

(2-15) ||Aw||L°°(GnQ r (o,o)) < ^ r °- 

Therefore, by using (12.151) in (12.141) we see that in G D Qi(0, 0) one has 
(2.16) \Av — v t \<C5. 

Furthermore, by (12.121) we obtain from (12.91) the following bound 

Iblll/xqGnQifO,!))) - 1 ‘ 

In addition, v vanishes on dG fl Qi(0,0). Therefore, if we let <5 —> 0 along a sequence, we will 
have by compactness (by using uniform interior H 2+a estimates and boundary H 1+a estimates) 
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that for a subsequence <5 —>• 0, we have v = v(5) —> vq uniformly on compact subsets, where the 
limit function Vo has the properties 


(2.17) 


A no — D t vo = 0 in 5^ x (—1,0], 

< ||vo||z,°° < 1, 

= 0 on ({x n = 0} n Bi(0)) x (-1,0]. 


Here, B+ = Hi fl {x n > 0}. We now denote by Vo the odd reflection in x n of the function 
vo to the whole -£>i(0) x (—1,0]. Then, A Vo — D t V o = 0 in F>i(0) x (—1,0], and Vo satisfies 
the remaining two conditions in (12.1711 above. In particular, from the smoothness of Vo up to 
({x n = 0} fl £?i(0)) x (—1,0] and the third property in (12.1711 above we see that Di Vo = 0 for 
i = 1,..., n — 1, and D t V o = 0 on ({ x n = 0} fl l?i(0)) x (—1, 0]. This gives Djj Vo = 0, Dj t \o = 0, 
i,j = 1, ...,n — 1, and D t tV o = 0 in ({ x n = 0} fl Hi(0)) x (—1,0]. Furthermore, since Vo is odd 
in x n we also have D nn V o(0,0) = 0. Using the fact that the variable t has weight two, from the 
Taylor expansion of Vo at (0,0) we obtain that there exists p > 0 universal (p < (VC) -1 ^ 1- ") 
would do), such that 

(2.18) 11vo — % n Qo \ (o) x [—p 2 ,o]) — Cp 3 < -^P + , 

where Qo(x) = E ”!] 1 qpxi + qo is an affine function which is approximating (note that q n = 0 is 
a consequence of D nn \ q( 0, 0) = 0), with 


(2.19) 


Ittl < C-i 


for some C\ > 0 universal. In particular, the product x n Qo(x ) is harmonic. We now fix such a 
universal p > 0. Since v = v(5) —> vq uniformly on compact sets as 5 —> 0 on a subsequence, by 
compactness we see that for 5 sufficiently small we have, 

1 

V 

From (12.2011 and (12.1811 we thus conclude that for 5 sufficiently small 

(2.21) ||v — £ n Qo|li/x,((5 n Q p ( 0i o)) < 11^ — v o|lioo( ( 5 n Q p (o i o)) + ll u o ~ x « < 3o|lx,<x>( ( 5nQ p (o,o)) 

1 1 


( 2 ' 2 °) V o|lioo(c n g p (o,0)) — 4 P 2+a - 


we have that 


+ ^ 


2+a 


\u x n\\ L o O ^G n Q p ^ !0 ^ < S. 


Moreover, from m and 

( 2 . 22 ) 

Thus, from the triangle inequality, if with C\ as in (|2.19l) we further restrict 6 in dependence of 
p so that nC\5 < \p 2+a , we have 

||v - uQ 0 \\ LOO (Q n Q p ( 0 ^ < p 2+Q . 


(2.23) 


The conclusion now follows by taking P(x) = P(x) + v 1+ "Qo(f) and by rewriting v, u in terms 
of v,u. We notice explicitly that (12.1011 follows from the definition of P and from the fact that 
(12.1911 gives 

||Qo(-)||L°°(GnQ r (o,o)) < n C- 

By using (12.1211 . (12.131) and the definition of P, we finally obtain 

11 v 'U'P\\ L oo (GnQp r (0,0)) = | |w 2+ "7)(—, - 2 ) - r 1+ “n(5o(-)||L°°(GnQ pr .(o,o)) 

= r 2+ “||v(-, -£) - «(-, ^)Qo(-)||L»(GnQ pr (o,o)) 

< ( pr) 2+a , 


r r* 


where in the last inequality we have used (|2.23[) . This proves (12. lip , thus completing the proof. 
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□ 


With this lemma in hand, one can establish Theorem 12.21 bv arguing as in [DSlj . We never¬ 
theless provide the details for the sake of completeness. 

Proof of Theorem \2.2l Now suppose there exists an affine function P for which the hypothesis 
of Lemma 12.61 holds for some r > 0. Then with v,u as in the proof of Lemma 12.61 we first note 
that by H l+a estimates up to the boundary as in Theorem 4.27 in [Li] and from the fact that 
v vanishes on G D Q 2 ( 0 , 0 ), we have for every (x, t) £ G D Q 1 / 2 ( 0 , 0 ) 


(2.24) 


|u(.x,t)| < C\d(x, G t ), 


where Gt = {x | (, x,t ) £ dG n <5i(0,0)}, and d(x,Gt) is the Euclidean distance of the point 
x £ M n from Gt- Note that (12.24p is a reformulation of the Lipschitz estimate in the x variable 
at the boundary. Moreover, from (12.51) and the definition (12.131) of u we find 

u v > c > 0 , 

where c is the universal constant in (12.51) . From this inequality we easily obtain 

(2.25) u(x,t ) > C 2 d(x,($t). 

The estimates (|2.24p and (12.251) imply that 

(2.26) \v\<Cu in GDQi / 2 (0,0). 

We now see how the desired conclusion (12.41) above can be derived from (I2.26|) and from 
Lemma 12.61 

First, by rewriting v,u in terms of v and u using (12.121) . (12.131) . we obtain as a direct conse¬ 
quence of (12.261) 

(2.27) \v - uP\ < Cur 1+a in G D Q r/ 2 (0,0). 

Moreover, since by (12.251) the function u is bounded away from zero from below in G(lQi/^(\e n , 0)), 
from H 1+a estimates for v and the fact that this function satisfies (12.161) . we also have 


(2.28) 


V 

I ~lltfi+«(GnQ 1 / 4 (§e n , 0 )) - C '(ll u llL oo (GnQ 2 (0,0)) + l ) ■ 


I v I 

' U "«I/4V2 l 

In (12.281) we have used the fact that, because of (12.81) and the definition of G, at each time level 
t £ (—1,0] the set K t = {x \ (x,t) £ G fl Q\//±(\e n , 0)} is at a Euclidean distance from Gt which 
is bounded below by C 5 , for some C 5 > 0 universal. In addition, the identity 


(2.29) 


u(x, t ) 

which follows from (12.121) . (|2.13|) . implies that 

(2.30) r " ,V 

u 

and 

(2.31) 


v ( x .*) =p( x ) + r l+a.^r’P) 


u(f,^y 

[D X (-)} H ^ G n Qr/4 (re n ,0)) = [Dx ( w)] H a (GnQ 1/4 ( , 0 )) 


V V 

< - >l+a;GnQ 7 . /4 (§e„,0)) = < - >1 


^ l+a;GnQi/4(2 e n ,0)) 

(See page 46 in [Li] for the definition of the seminorm < w >i+ Q; n which corresponds to the 
i^p-Holder seminorm of w in t). 

Summing up, (12.281) . (I2.29|) . (12.301) and (12.311) imply the following estimate 
IHItf 1+ “(GnQ r/4 (§en,0)) < C'(ll^llL oo (GnQ 2 (0,0)) + l) • 


(2.32) 
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We would like to mention that the region GrK5 r / 4 (|e n , 0) is at a parabolic distance proportional 
to r from dG, a fact that follows from the first inequality in (12.81) which represents a flatness 
assumption on the boundary. 

Multiplying v by a suitable constant, we may now assume that the hypothesis of the Lemma 
12.61 is satisfied for ro small and P = 0. Given any 0 < r < 1, with p fixed as in Lemma 12.61 we 
choose k G N such that p k+1 r q <r< p k r o- We thus apply Lemma 12.61 iteratively, i.e., first for 
ro, then for pro, p 2 ro and so on. We finally obtain a limiting affine function Po such that 

(2-33) ||i> - uPq ||L°°(GnQ r (o,o)) < Cr 2+a for r < r 0 . 


We note explicitly that 

OO 

p o(x) = ^2(p l ~ lr o) 1+a Qi(^]—), 

P r o 

where Qi is the affine function obtained after the i-th application of Lemma 12.61 in the iteration 
argument. Given that (12.331) holds with P = Po for every r < ro, therefore for any given r < ro, 
we can now repeat the arguments which lead to (12.271) with P = Po and consequently obtain for 
all (x,t) G GnQi(0,0) 


(2.34) 


u 


-(x,t) — Po(x) <C(\x • 


l+a 

) 2 . 


This implies the iL 1+Q -regularity at the boundary point (0,0). Combining ()2.34l) with (12.32[) . 
the desired conclusion follows by arguing for instance as in the proof of Proposition 4.13 in |CC] , 
We mention that, although the latter result is for the elliptic setting, the same argument goes 
through in the parabolic setting when the Euclidean distance in M ra is replaced by the parabolic 
distance in M n+1 . 

□ 


3. Higher regularity 

In this section we establish the case k > 1 of Theorem E2J Our main result is the following. 

Theorem 3.1. Let G be of class H k+a in G fl <52(0,0), the other assumptions on G being as 
in Section 0 Let u and v be two solutions of the heat equation in G D <52(0,0) such that u,v 
vanish on d p G D <52(0,0). Also, suppose that u > 0 in G 0 <52(0,0) and assume that it satisfy 
the normalization m- Then, for some C > 0 universal one has 

V 

(3-1) ll“ll// fc +“(GnQi(o,o)) ^ C (11^1 lz,°°(G'nQ 2 (o,o)) + l) • 

Before proving the theorem above, we introduce some additional notations. Henceforth, we 
let No = N U {0}. Given a multi-index m = {m\, G Nq, for x = (. x\,...,x n ) G W 1 we 

denote by x m = x™ 1 ... x r ff n the monomial of degree \m\ = mi + ... + m n . Henceforth, for 
i = 1, ...,n, the notation i will indicate the multi-index in Nq which has 1 in the i-th position 
and 0 elsewhere. 


Definition 3.2. By a parabolic polynomial of degree k G No we mean an expression of the form 

P{x,t) = a m,ex m t e , 

0<\m\+2e<k 

where a rn j G M. Given such a P we define its norm as follows 

||P|| = max | a m A. 

0<|m|+2t<fc 

Following jDSlj. see also Chapter 1 in m, we introduce the following dehnition. 
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Definition 3.3. We say that a function f is pointwise H k+a at (0,0) if there exists a parabolic 
polynomial P(x,t) of degree k such that f(x,t ) = P{x,t) + 0 ^(|x[ J + . We indicate this 

with f £ H k+a (0,0) and we denote with ||/||jyfc+cqo,o) the smallest M > 0 such that ||P|| < M 

and | f(x, t) — P(x, f)| < M{\x\ 2 + for t < 0. We say that f is pointwise H k+a at (xo,to) 

if the function h(x,t) = f(x + xo,t + to) is pointwise H k+a at (0, 0), and we indicate this with 
f £ H k+a (xo,to). Finally, given a bounded set G C M n+1 we say that f £ H k+a (G) if 

de f 

\\f\\H k + a (G) = SU P \\f\\H k + a {xo,t 0 ) < °°- 
(■ xo,to)£G 

We note that the class H k+a (G ) coincides with that defined on p. 46 in [Li] , 

Now with u, v as in Theorem I3J1 it follows from the Schauder theory (see Chapters 4 and 5 
in M) u -> v are i n H k+a up to G n < 32 ( 0 , 0). As in Section [2 we assume as well that 

(3.2) p = p 0 (0) = e n , f (0,0) = 0, D'f( 0,0) = 0, and \\f\\ H k +a < c 0 , 

where cq is a dimensional constant which is chosen sufficiently small in such a way that (e n , —3/2) £ 
G. This latter property can always be achieved by suitably scaling the domain as in (I2.6D - (12.8|) . 
Furthermore, we assume that 

(3.3) Du{ 0,0) =e n , \\u - x n \\ H k+ a ^ GnQ2(m) < 6, ||/||^+a<A 

The second inequality in (13.31) can be seen as follows. Since u( 0,0) = 0 there is a parabolic 
polynomial P u , of degree at most k, such that for all (x,t) £ G D Q3/2(0,0) one has 

(3.4) \u(x,t) - P u (x,t) | < C(|x| 2 + l^l) - ^. 

Note that such a polynomial corresponds to the weighted Taylor expansion of order k of u at 
(0,0) in which derivatives in t are assigned the weighted order 2. If we let u rQ be as in (12.61) . 
then with Du{ 0,0) = e n we have that 

r 2 

(3.5) \u ro (x,t) — x n \ < —\Pi(r 0 x,rQt)\ = r 0 \Pi{r 0 x, r%t)\, 

ro 

where P\(x, t ) = P u {x, t) — x n . Therefore, by choosing ro sufficiently small the second inequality 
in (13.31) can be ensured with our new u = u ro and G = G r ° , following computations similar to 
those in Section [2j 

Let P be a given parabolic polynomial of degree k. By using the fact that A u — ut = 0, we 
have that 


(3.6) (A — Dt)(uP) = 2 < Du, DP > +u(AP — P t ). 

As a first step we write a formula for (13.61) when P(x,t) = x m t l , for a given multi-index rri £ Nq 
and t £ No- This is (13.81) below. To derive such formula we first notice that from (13.31) it follows 
that u = x n + w where w is of parabolic order > 2 and ||«>|l.ff fc -i-“(G n Q 2 (o, 0 )) — O n the other 
hand, (13.41) gives u(x,t) = P u (x,t) + z(x,t), where P u is a polynomial of degree at most k. 
Combining these two facts, and letting P\(x,t) = P u {x,t ) — x n , we can thus write 

(3.7) u(x,t) = x n + P\{x,t) + z(x,t), 

where, we note explicitly, the polynomial Pi is of degree at least two, provided it is nonzero. 
We thus have 

Du(x, t) = e n + DPi(x, t ) + Dz(x, t). 

This gives 


2 < Du, DP > +u(AP — P t ) = 2 D n P + x n (AP — P t ) + 2 < DP \, DP > 

+ 2 <Dz,DP> +(Pi +z){AP-P t ). 
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Keeping in mind that 

D n P = m n xx n D nn P = m n (m n - 1 )x m ~ n t l , 


we find 

(A - D t )(uP) = m n (m n + 1 + Y - l)x m ~ 2UjT t e - 

i^n 

+ {2 < DPi, DP > +Pi(AP - P t )} + {2 < Dz, DP > +z(AP - P t )}. 

It is now easy to see that 

2 < DP\, DP > +P 1 (A P -P t )= Y c^x q t K . 

0<M+2K<|m|+2£+fc-2 

We note explicitly that, since P\ is of degree at least 2, one has that c ™’ K f ^ 0 only when 
\m\ + 21 < |g| + 2 k. However, for later purposes we find it convenient to isolate from the right- 
hand side of the latter equation a parabolic polynomial which is of degree at most k — 1, i.e., we 
decompose 


E 

0<|g|+2K<|m|+2£+fc-2 


u q,K ^ L 


E 

\m\-\-2i<\q\+2n<k —1 

+ ~E~ 

k<\q\+2n<\m\+2t+k-2 




m,l q + K 
u q,K ^ 1 ' 


If we define 


w m j{x, t) = 2 < Dz, DP > +z(AP - P t ) + Y c ™K X<lfK - 

k<\q\+2K<\m\+2£-\-k—2 

then we conclude that 

(3.8) (A - D t )(uP) = m n (m n + 1 )x m ~ Yi t i + Y ™i(™i “ 1 )x m ~ 2i+Ji t l - 

i^n 

+ Y c 7^x q t K + w m j(x,t), 

\m\+2£<\q\+2K,<k—l 

where we have | Cq'£ | < C8, a fact which follows from | |i.o, (c?nQ 2 (o,o)) — Moreover, again 

from (13.31) we have for all (,x , t) G G fl Qi(0,0) that 

(3.9) \w m /(x,t)\ < C5((\x\ 2 + |t|) 2 + ). 

In (13.91) we have crucially used the fact that \Dz(x,t)\ < C<5(|.t| 2 + \t\) 2 ^) for all (x,t) £ 

GnQi(0,0). 

We now turn to determining a suitable expression for the right-hand side in (13.61) in the case 
in which 

P(x,t) = Y a m/x m t t 

0<\m\+2i<k 

is a general parabolic polynomial of degree k. In such case, we obtain from (13.81) 

(3.10) (A - D t )(uP) = R{x)+ w(x), 
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where R is the parabolic polynomial of degree k — 1 given by 

(3.11) 

R(x,t)= ^2 dq jK X q t K = ^2 a m,e 

0<|g|+2re<fc—1 0<\m\+2t<k 

- ix m+ji t 1 - 1 + ^2 c ^ xH 

\m\+2£<\q\+2K,<k—l 

and 

w(x) = ^2 a m,eWmA X )- 

0<\m\+2£<k 

Note that (J3.11D gives 

(3.12) ( q n + 1 ){Qn + A a q+n,K + ^ T 1 ){Qi + ‘2,) a q+2i—n,K — ( K + l)o<j— n,K+l + C q ’ K a m y = 

i^n 

Given a parabolic polynomial of degree k — 1 such as i2(x) = X)o<M+ 2 K<fc-i d q , K x q t K , our ob¬ 
jective is finding a parabolic polynomial of degree k, P(x,t) = X)o<|m|+ 2 f<fc a m,e xm t, such that 
(13.101) hold (in particular, we will be interested in this section in the case when R = 0, see Defi¬ 
nition [33] below). This will be possible if, given constants d q>K , we can solve the linear system 
(13.121) above for the unknowns a m (_. Notice that one can think of (13.121) as an equation where 
a q+ n tK is a linear combination of d 9)K ’s and a m /s such that either \m\ + 2£ < |g| + 2k + 1, or when 
\m\ + 2l = q + 2k + 1, then m n < q n + 1. It thus follows that we can solve (13.121) if we arbitrarily 
assign all the coefficients a m j when m = (mi,... ,m n ) is a multi-index having m n = 0. In this 
respect we emphasize that the crucial fact which makes this claim possible is that the third term 
in the left-hand side of (13.121) . namely, the one which comes from differentiating in the time 
variable t, has the same degree as the first two terms. We define the order of a coefficient a mt g 
as \m\ + 2£, i.e., the weighted degree of the corresponding monomial x m t e . 

To verify the above claim we briefly describe the procedure of determining the coefficients. 
First of all, O(o,...,o),o> which is the unique coefficient of order 0, is assigned arbitrarily since it 
trivially satisfies the requirement m n = 0. We proceed by a double induction. Suppose we know 
all coefficients a m j up to order p. Given a coefficient a mo ^ 0 of order p + 1, i.e., |mo| + 2^o = P + 1, 
let (mo)n denote the entry at the n-th position of the corresponding multi-index (mo, £q)- Clearly, 
0 < (m 0 )n < p + I- We determine all coefficients of order p + 1 by induction on (mo) n ■ First, 
all coefficients a mo ,£ 0 of order p + 1 with (mo) n = 0 are arbitrarily assigned. Suppose now all 
coefficients a mo j Q of order p + 1 with (mo) n < x are known. Then, given a coefficient a mi ^ of 
order p+1 with {m\) n = y +1, we have from (|3.12l) that such a coefficient is expressible in terms 
of lower order coefficients a m j. which are already determined, and coefficients a mo ^ 0 of order 
p+1 such that (wio) n < x an d which are supposed to be known by the induction hypothesis on 
(m 0 )n- Therefore, all coefficients a rnQ of order p + 1 with (mo) n = x + 1 can be determined 
once all coefficients a mo j 0 up to order p+1 with (mo) n < x are known. In this manner, all 
coefficients up to order p+1 can be determined once all coefficients of up to order p are known. 
The claim thus follows. 

For a fixed k as in Theorem 13.11 we now consider parabolic polynomials of degree < k. We 
next introduce a notion which generalizes to the case k > 2 that in Definition 12.51 above. 

Definition 3.4. We say that P(x,t) = X)o<|m|+ 2 ^<fc a m,l x ' m "P an approximating parabolic 
polynomial of order k for ^ at (0,0) if we have R(x, t) = 0 in the representation (13.101) above. 
This is equivalent to the fact that the a m y satisfy (13.121) with d q , K = 0. 

Remark 3.5. The motivation for Definition \3.f \ comes form the fact that, when P is approx¬ 
imating, then from (13.101) we obtain H{uP) = w, with w of weighted order k — 1 + a. This is 


m n (m n + l)x m n t e 


- l)x m ~ 2l+n t e 


i^n 
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crucially used in (13.1811 below since it allows us to cancel the term r k 1+a from both sides of the 

_l - 1 - Q! 

equation. We have in fact, see (ESI), w(x,t ) = (|x| 2 + \t\) 2 w\(x,t), where ||w;i(x, t)\\ < 5. 

Note that when k = 1, and therefore P(x, t ) = ao + X^=i a i x i> this notion is equivalent to saying 
that x n P is caloric, which is in turn equivalent to the condition a n = 0. In that case from 
(12.151) we see H(uP ) < C5r a , which in the case k = 1 is the w of order 1 — 1 + a = a. We 
also mention that in the case of the variable coefficient operators treated in Section [7] we will 
no longer impose the condition R = 0 in the definition of approximating polynomial. This is so 
because of a nonzero right-hand side. 

The proof of Theorem 13.11 follows the same lines as that of Theorem 12.21 once the following 
lemma is established. 


Lemma 3.6. Assume that for some r < 1 and P an approximating polynomial of order k for ^ 
at (0,0) with ||P|| < 1, one has 

(3-13) ||v - 'wP|| L oc (GnQr . ( o i o)) < r k+1+a . 

Then, there exists an approximating polynomial P of order k such that for some C, p > 0 
universal, we have 

(3-14) \\P-P\\L~(Gn Qr m)<Cr k + a , 

and 

(3-15) ||u - uP\\ L ^ {GnQpr{m) < ( pr) k+l+a . 


Proof. With G as in the proof of Lemma 12.61 we define v(x,t) by the equation 


(3.16) 
and we let 

(3.17) 


v(x,t) = u(x,t)P(x,t ) + r k+1+a v(~, -A), 


x t 


r r" 


u(x,t ) = 


u(rx, r 2 t ) 


Since v is a solution of the heat equation we have 


x t 


0 = Av - v t = (A - D t )(uP) + r k 1+a (A7 - D t v)(~, -^). 




Using the fact that P is an approximating polynomial of order k and (13.101) . we conclude 


(3.18) r k ~ 1+a (Av - D t v)(^ ±) = -w(x, t). 

We emphasize the crucial role played in the latter equation by the property of P being an 
approximating polynomial, see Remark 13.51 above. By (13.91) we conclude that in G D Qi(0,0) 
one has 

(A - D t )v = h, 

where 

h e H k ~ 2+a (G D Qi(0,0)) and \h\ < CS. 

In addition, v vanishes on dG D Qi(0,0). Therefore, if we let 5 —> 0 along a sequence, we 
will have by compactness (by using uniform interior H k+1+a estimates and boundary H k+a 
estimates) that for a subsequence 6 —> 0, we have v = v(6) —> vo uniformly on compact subsets, 
where the limit function vg has the properties 


Avq — D t v o = 0 in B+ x (—1, 0], 

( | |^o | |l°° A 1) 

^v 0 = 0 on ({x n = 0} D Bi(0)) x (-1,0]. 


(3.19) 











A PARABOLIC ANALOGUE OF THE HIGHER-ORDER COMPARISON THEOREM, ETC. 


13 


Proceeding as in the proof of Lemma 12.61 we can find a polynomial Q o of degree k such that 
x n Q o solves the heat equation, and 

(3-20) lb - ZnQolL<x, (s + x( — P 2 i0]) < Cp k+2 < \p k+1+a , 

for some C, p > 0 universal. This follows for instance by odd reflection of vo across x n = 0 
(which provides again a solution to the heat equation), and by applying Lemma 1.1 in m to 
the extended function. Then, again for 5 > 0 sufficiently small, we have 

(3- 21 ) 11^ ~ ^ollicxj^nQpfo.o)) — ~^P k+1+a - 

Therefore, from (13.201) . (13.211) and the triangle inequality we obtain, 

(3.22) \\^ ~ x nQo\\L°°(GnQ p (o,o)) — 2^ ' ++ ' 

Now, by using the second inequality in (13.31) we conclude that, if the choice of 5 is further 
restricted in such a way that CS < \p k+1+ot , for some C which depends on ||Qo|| and hence is 
universal, then 

Hu — «Qo|lioo((5nQp(o,o)) — 4 p k+l+a ■ 

Therefore, by rewriting v and u in terms of u and v we have that 

(3.23) ||„ - „(p + r‘+«Q„(-, ^))ll L «(o nQ , r( o, 0) ) < jW t+1+ “. 

At this point, (13.231) would complete the proof of the lemma if we knew that the parabolic 
polynomial P + r k+a Qo (^, - 7 ) is approximating for Unfortunately, this is not necessarily 
the case and we need to further modify Q 0 to some other polynomial Q, without essentially 
modifying the estimate (13.231) . 

Let us write Q 0 (x,t) = Eo<M+2«<k b q , K x q t K and Q(x,t) = Eo<| 9 |+ 2 «<fc \ K x q t K . In view of 
(13.121) . and by using the fact that P(x,t ) = Eo<|m|+ 2 £<A; a m,ex m t £ is an approximating polyno¬ 
mial for 0 of degree k, we obtain that, in order for 

(3.24) P(x,t)= f P(x,t)+r k+a Q (- 1) 

to be an approximating polynomial for ^ of degree k, the coefficients b QyK of Q should satisfy 

(3.25) (q n + 1 ){q n T 2')bqj r n )K + ^ ^ (qi + 1 ){Qi T ‘^‘)^ ) q+ 2 i—n,K — ( K T 1)^ 9 —n,re+l + c q : k bm,£ = 0, 

i^n 

where 

pn,t \q\+2K+l-(\m\+2£) m,£ 

The justification for (13.251) is as follows. If we write P(x, t)+r k+a Q( 4j-) = Eo<| 9 |+ 2 «;<fe ^q,K.x q t K , 
then it is clear that 

(3.26) a qiK = a q , K + r fc+Q_|g|_ 2 K 6 5A 

Now, imposing that P(x, t) + r k+a Q(f, \) be approximating of order k for ^ at (0,0) is, in view 
of (13.121) and Definition 13.41 equivalent to 

(3.27) ( q n + 1 )(<?n 4" ^ ' (<& T 1 )(.Qi T 2)hq_|_2i_ji iK ~ (k + l)o 9 — n,K+l + Cg ’ K = 0. 
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Replacing (13.2611 in (13. 2711 . and using the fact that P is itself approximating for ^ (see (13.1211 1. 
after some elementary computations we recognize that the coefficients b ,must satisfy 


„k+a 


r -(\ q \+2 K+ i) {qn + 1)(9n + 2 )b q+ n, K + r-d 9 l +2K+1) + 1 )(9i + 2)6 9+2? _ nA 

i^n 


- r-M+ 2K+1 \K + l)6 g _n,K+i + r-CI^I+^c^^ 


= 0. 


Eliminating r k+a and then multiplying the resulting equation by r kl+ 2K +i we finally obtain 
(I3.25H . Recalling that 0 < r < 1, and that from the defintion of d Q)K in (I3.11|) one has c^f / 0 
only when \m\ + 21 < |g| + 2k < k — 1 , we conclude that 

|c^| = r ,kl+2«+l-(|m|+2Q| c m^| < r | c m^| < ^ 

Since (13.2011 above implies that Q o is approximating for j 2 -, then its coefficients b m ^ must satisfy 
(3.28) (q n + 1 )(<7 n + 2)b q +n,K + ^^(<7i + 1) + 2 ) 6 g+ 2 j_ 7 j ]K — (k + l)&g_n jK +i = 0. 

i^n 


If we now subtract (13.2811 from (13.2511 , we find that the coefficients of Q — Q o solve a linear system 
with left-hand side bounded by C5 and contains unknown coefficients b rn £ such that \m\ + 2£ is 
less than the sum of the indices of the coefficients of Q — Q o appearing in the right-hand side. 
We in fact have 


(3.29) —c^j^bmj — (■ q n + l)(q n + 2)(b q+ n,K — b q+ n )K ) 

+ + !)(?* + 2)(bq+2i-n,K, ~ ^q+2i-n,K,) — ( K + l)(&q-ra,K+l — ^q-n,K+ 1)- 

i^n 

The reader should note that the order of the coefficient of any term in the right hand side 
of the latter equation is |g| + 2k + 1 > \m\ + 21. Again, this is so since c q .n / 0 precisely 
when \m\ + 21 < |q| + 2k. Consequently, if we set b m j = b m j when m n = 0, then in view of 
the procedure described after (13.12|) . we can determine all the other coefficients b q ^ K of Q by 
induction on the order of the coefficient \q\ + 2k. Morever, since the coefficients of Q — Qo solve 
a linear system with left-hand side bounded by C5 , we can further ensure that 


IIQ “ Qo||l°°(Qi(o,o)) < C5. 

Since 

llQll.L°°(GnQp(0,0)) — I IQ “ Qo|li,°°(GnQp(0,0)) + IIQIlL°°(GnQp(0,0))’ 
we conclude that there exists a universal constant C > 0 such that 


HQIlL°°(GnQ p (o,o)) — C- 

Therefore, by (|3.22l) and by choosing a smaller 5 if needed, one can ensure that 

(3.30) ||h — x nQ|| i oo( ( 5 nQ p (o,o)) — II® ~ ^nQo||£oo(c n Q p( - 0j0 )) + \\x n (Q — Qo)||£°o(G n Q p ( 00 )) 

< \p k+1+a + C5p < | p k+1+a . 

Then, again by the second inequality in (13.311 and by (13.301) . we obtain for a smaller choice of S 
if needed that 


(3.31) ||h '“QIIlo^GiXMo.o)) — 11® Xri QllL°°(GnQp(o,o)) + IIQ(® x «)llL°°(GnQp(o,o)) 

< 5/+i+« + < 2/ + ‘+" + cs < / +1+ “ 
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Recalling the definitions (13.161) . (13.171) and (13.241) . the conclusion (13.151) of the lemma is now 
obtained as follows 


\\v ~ uP\\ L <*>( G nQ pr m) = II uP + r + +a v(~^) ~ <P + r +a uQ(-, ^)\\ L ~ (GnQprm) 

= r k+ 1 +a \\v(-, - u(-, ^)\\L°°(GnQ pr m) 

= r k+1+a \\v - uQ\\ L ~ {G n Qp m) ^ (. pr) k+1+a , 
where in the last inequality we have used (13.311) . This completes the proof of the lemma. 


□ 


Proof of Theorem, \3 . 1\ The proof of the theorem now follows by iterating Lemma 13.61 To start 
the process of iteration, we take P = 0. Multiplying v by a suitable constant, one can ensure 
that the hypothesis of Lemma 13.61 holds for some ro universal. The rest of the proof remains 
the same as in the case k = 1, in which Lemma 13.61 is applied iteratively first for ro, then for 
pro, p 2 ro, and so on. We finally obtain a limiting polynomial Po °f degree at most k having the 
following representation 


Po{x,t) = Y J (p i ~ 1 ro) k+a Qi( 

i= 1 


X t 

p^ro' (p* -1 ?’o) 2 


where Qi is the polynomial obtained after the f-th application of Lemma 13.61 Furthermore, the 
following holds 

(3-32) \\v -uP 0 \\ L °°(GnQr(.o,o)) <Cr k+1+a , r<r 0 . 


For r > 0 we now consider the functions v,u defined, with P = Po, as in the proof of Lemma 
13.61 Given that (13.321) holds for any given r < ro, we can argue as in the case k = 1 and obtain 
that 


(3.33) \v\<Cu in G fl Q 1/2 (0,0). 

By rewriting u,v in terms of u and v, from (13.331) we obtain for (x,t) £ G D Q r /2(0,0) 

(3.34) | v(x,t) — u(x,t)Po(x,t)\ < Cu(x,t)r k+a . 


The inequality (13.341) gives for (x,t) £ G D Qi(0, 0) 


I ~{x,t) - Po(x,t )| < C(|x| 2 + |7|) fc 2 “ - 
u 

This implies P fc+Q -regularity at the boundary point (0,0). Finally, by arguing as in the case 
k = 1, the identity 


(3.35) 


h^ = F„(x, i)+ 

u{x,t) u(f, 


1 V 
71) 


implies that the H k+a norm of ^ is bounded in regions of the form G D Q r / 4 (|e n ,0), i.e., in 
regions which are away from the boundary of G by a parabolic distance proportional to r. At 
this point, the conclusion follows similarly to the case k = 1 by arguing as in the proof of 
Proposition 4.13 in |CC| . 


□ 
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4. VARIABLE COEFFICIENTS 

In this section we intend to extend Theorem 12.21 and Theorem 13.11 to variable coefficient 
parabolic operators of the type 

n n 

(4.1) Lu = CLjj(x, t)DjjU + h(x, t)DjU + c(x, t)u — ut = 0, 

*,.7=1 »=i 

satisfying appropriate regularity assumptions on the coefficients, and to strong solutions u £ 
Wp’i oc (G). Henceforth, we use the notation Tr(M) for the trace of a matrix M. We can thus 
write Y^ij=i a ij{ x i t)Dij u = Tr (A(x,t)D 2 u), where we have indicated with D 2 u = [D t ju] the 
Hessian matrix of u. 

As in the case of the heat equation, in order to better present the ideas we first treat the case 
k = 1 in Section RTTl Subsequently, in Section l4~2l we treat the case of k > 1. In the sequel we 
will denote with Wp’ 1 (G) the parabolic Sobolev spaces. For their precise definition we refer the 
reader to p.155 in [Li] . We will indicate with W^f oc (G) the standard local spaces. 


4.1. H 1+a regularity. The assumptions on G are as in the hypothesis of Theorem 12.21 In this 
section we assume that A = [a,j] £ H a (G), b,c £ L°°(G). The following is our main result. 

2 1 

Theorem 4.1. Let p > 1 and suppose that u £ W p ’ loc (GC\ <^ 2 (0, 0)) be a positive strong solution 
to (14.11) above. Let v £ W P ' 1 {G D Q2(0,0)) be a strong solution to 

(4.2) Lv = g , in G , 

where g £ H a (G D Q 2 (0, 0)). Assume that u,v vanish on d p G D Q2(0,0). Furthermore, let u 
satisfy the normalization condition (USD. Then, one has 

V 

(4-3) || -||lL 1 +“(GnQi(0,0)) < C'dl' (; llL° o (GnQ 2 ( 0 , 0 )) + IMI//“(GnQ2(0,0)) + 1). 

for some C > 0 universal. 


Remark 4.2. We first note that from the assumptions on the coefficients, the Calderon- Zyg- 
mund theory implies that u £ W^’^ oc (G) for all 1 < q < 00, see for instance Proposition 7.If in 

P. Then, we can invoke Theorem f.29 in E3 to conclude thatv,u are in H l+a {G fl < 32 ( 0 ,0)). 
We note that Theorem f.29 in [Li] can be applied to strong solutions in loc via approxima¬ 

tions by solutions to equations with smooth coefficients and by an application of the comparison 
principle Theorem 7.1 in [Li]. We refer to G0 for the elliptic counterpart of such intermediate 
Schauder type regularity result. 


After a suitable change of coordinates and parabolic dilation similar to (12.61) . we may assume 
that 


(4.4) 


'A(0,0) =/, 

/ (0) = 0, T>7(0) = 0 , ||/|| C 1. a (GnQ2 (0,0)) — ^5 

Du{^ 0,0) = e n , | \u — x n \ |// 1 +Q!(GnQ 2 (o,o)) — ^ 

max < [A] aiGn Q 2 ( 0i o), lbllif“(GnQ 2 (o,o))dl^i c llL°°(GnQ 2 (o,o)) \ < 


Here, [A] a indicates the a-Holder seminorm of [ 0 ^], see p. 46 in [Li], More precisely, first by a 
suitable change of coordinates, we can ensure that A(0, 0) = /. Then, by letting 


u(r 0 x, rgt) 
ro 


v(r 0 x, rgf) 
ro 


(4.5) 


u ro (x,t) 


v ro (x,t) 
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as in (12.611 , we have that u ro , v ro solve in G r ° (same definition as in Section [2]) 


where 

and 


Lqu To — 0 , Lqv To — go, 

Low = Tr(AoD 2 u;)+ < bo, Dw > +cqw — wt, 


Ao{x,t) = A(r 0 x,rlt), b 0 {x,t ) = r 0 b(r 0 x, r%t), c 0 = roc(r 0 x, r^t), g 0 (x,t) = r 0 g(r 0 x,r%t). 

Therefore, if ro is suitably chosen depending on S, (14.41) can be ensured. As before, by abuse 
of notation, we keep calling u ro = u, G r ° = G and so on. Moreover, as in Section [21 5 will be 
determined later. 

We now introduce the relevant notion of approximating function with respect to L and g, 
where g is as in Theorem 14. 11 


Definition 4.3. We say that P(x) = ao + aiXi an approximating affine function at 
(0,0) for ^ with respect to L,g if 2a n = <7(0,0). 

With this notion the corresponding statement of Lemma 12.61 remains the same, but its proof 
needs to be slightly modified. 


Lemma 4.4. Let u,v be as in Theorem \f.l\ Assume that for some r < 1 and P(x) = oo + 
££=1 a i x i an approximating affine function at (0,0) for (with respect to L and g) with \ai\ < 1, 
one has 

(4-6) \\v - uP\\ L °°(GnQ r (o,o)) < r 2+a . 

Then, there exists an approximating affine function P such that for some C, p > 0 universal, we 
have 

(4-7) \\P ~ P\\L~{GnQ r m) <Cr 1+a , 

and 

(4-8) \\v - mP||z/x> (GnQpr(0i0)) < (pr) 2+a . 

Proof. We point out the essential modifications in the proof of Lemma [2. 6 1 in the present context. 
Let v and u be as in the proof of Lemma 12.61 Then, one has 

(4.9) g = Lv = L(uP) + r a Lv(x/r,t/r 2 ), 
where 

(4.10) Lv = Tr(AD 2 v) + r <b, Dv > +r 2 cv — vt, 
and 

(4.11) A = A(rx,r 2 t), b = b(rx,r 2 t), c = c(rx,r 2 t), ( x,t)€G. 

Since Lu = 0, one has 

L{uP) = 2 < ADu, DP > +u <b, DP > . 

From (|4.4I) we have for all (x, t) G GflQ r (0,0) 

' A(x, t) = I + 5r a M(x, t), 

Du(x,t) = e n +5r a w(x,t), 
g(x,t) = g(0,0) + 6r a wi{x,t), 
k | u(x, t)| < C5r, 

where M, w, w\ are bounded. Therefore, using ()4.12l) and (14.41) and the fact that P is approxi¬ 
mating for -, we obtain for some bounded function I\{x,t), 

\L(uP) — g\ = |2 < ADu, DP > +u<b,DP> -g\ = \g(0, 0) - g(x, t) + 5r a K(x, t)| < Kffiffi. 
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Combined with (14.91) this estimate gives 

\Lv\ < C6 in GnQi(0,0). 
We also note that, with u as in (12.131) . we have 

Lu = 0. 


Letting 5 —> 0, from (I4.4jl we obtain that up to a subsequence 5 —> 0, v = v(5) —> vo, which solves 
(12.171) . Note that, unlike the case of the heat equation, we do not presently have uniform interior 
H 2+a estimates for v. Nevertheless, because of H 1+a estimates for v = v(5) up to d p Gr\Qi( 0, 0) 
independent of 5, by applying Theorem 6.1 in jCKSj we can ensure that vq is a L p viscosity 
solution of the heat equation in the sense of [ CKS] . The regularity theory for viscosity solutions 
now ensures that vq is a classical solution of the heat equation. As in the proof of Lemma 12.61 
we have that (I2.18l) - (|2.23|) holds for Qo(x) = Y17=i ( L X % + Qo, with q n = 0. Then, as in the case 
of heat equation the conclusion of the lemma follows with P = P + r 1+a Qo{f)- Note that, if we 
let P(x) = Y17=i a i x i + a o and P{x) = X^r=i ®i x i + then since P is an approximating affine 
function for L and g we have a n = . Since q n = 0, we have 


a n. — T f Qn 


g{ o,Q) 
2 


This shows that also P is an approximating affine function at (0, 0) for - with respect to L and 
g. From this fact, the verification of (14.271) above follows as that of (12.111) in Lemma 12.61 


□ 


Proof of Theorem We repeatedly apply Lemma 14.41 To start the process we first take 
P(x) = ^Y^-x n . Then, (14.61) holds for some universal r = ro when v,g and P are multiplied 
by suitable constants. As before, by iterating Lemma 14.41 with r = r$, pro, p 2 ro and so on, we 
obtain a limiting affine function Po such that (12.331) holds. We note that in this case, Po has the 
following explicit representation 

OO 

(4.13) P 0 (x) = P(x) + 

2=1 

where P(x) = x n , and Qi{x) is the affine function obtained in the i-tli iteration of Lemma 
1441 The rest of the proof remains the same as that for the heat equation. 

□ 

4.2. H k+a regularity for k > 2. In what follows the assumptions on G are as in Section [3j 
We have the following higher-regularity result for variable coefficient operators. We note that 
in the next result we do not assume that u and v are strong solutions as in Theorem 4.1 since 
by the regularity theory one infers that both u and v are classical solutions. 

Theorem 4.5. Let u and v be (classical) solutions in G D Q2(0, 0) of the equations 

(4.14) Lu = Tr(AD 2 u)+ < b, Du > +cu — ut = 0, 
and 

(4.15) Lv = g, 

where A,ge H k ~ 1+a (G fl < 32 ( 0 ,0)), b, c £ H k ~ 2+a (G PI < 32 ( 0 ,0)). Assume that u,v vanish on 
d p G D < 32 ( 0 , 0 ). Also, let u > 0 in G D <32(0,0), and assume furthermore that it satisfy the 
normalization condition m- Then, one has 
v 

(4-16) ||-||^ + «(GnQi(o,o)) — C(lrllL oc (GnQ 2 (o,o)) + lblli?fc-i+“(GnQ 2 (o,o)) + 

for some C > 0 universal. 
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As before, from the Schauder theory as in Chapters 4 and 5 in U, we have that u,v G 
H k+a (G PI Q 3 / 2 (0,0)). By a suitable change of coordinates and parabolic dilations similar to 
(HU), we can assume that 


A(0, 0) = /, ||A — 7||.tf fc - 1 +“(GnQ 2 (o,o)) — ^ 

/(0) = 0, 147(0) = 0, ||/||H fe +“(GnQ 2 (o,o)) — ^ 

Du{ 0,0) = e n , \\u - x n \\ H k+ a{GnQ2{0fi)) < 5, 

max < ||5 r llH fc - 1 +“(GnQ 2 (o,o))) ll^> c llH fe - 2 +“((Gno 2 (o,o)) \ ^ ^ 


(4.17) 


where 0 < 6 < 1 is to be chosen appropriately later. 

Similarly to what was done in the proof of Theorem 13.11 above, we need to compute L(uP ), 
where P(x, t) = ]Co<|m|+ 2 £<fc a m/ xm t £ is a parabolic polynomial of degree k. Since by (|4.14l) we 
have Lu = 0, we obtain 

(4.18) L(uP ) = uLP + 2 < ADu, DP > + < b, DP > . 

By the linearity of L we are thus led to understand (14.181) when P(x,t) = x m t £ . In such case, 
from (14.171) it follows that 

(4.19) L[uP) = m n (m n + m^rm - l)x m “ 2 ~ i+Ji t e - tx m+ "t e - 1 

i^n 

+ Y, c™ k £ xH« + w m/ , 

\m\+2£<\q\+2K<k-l 

where because of the representation (13.41) which is valid for u, (14.171) and (14.181) . we have that 

(4.20) \<%£\ < CS, \w m/ \ < C5{{\x\ 2 + and \\'WmA\ H k - 2 +°‘(Gn Qr ) < C5r - 

In (14.191) . we used the fact that although the term u < b,DP >G H k ~ 2+a (G D Q 2 (0,0)), but 
nevertheless (14.191) and (14.201) can be justified as follows. We first note (see also (13.71) above), 
that because of (14.171) . we have 

(4.21) u(x,t) = x n + Pi(x,t) + wi(x,t), 
where P\ is a polynomial such that 2 < deg(A) < k and 

k-\-oc 


l-Pill < C5, |u;i(x,t)| < C6{\x r + \t\y a 


We furthermore note that, since b G H k 2+a (G D 72 ( 0 ,0)), we can write 

(4.22) b(x,t) = Pb(x,t) + bi(x,t), 

where Pi, is a vector field in M n each of whose components are polynomials of degree at most 
k — 2. Moreover, because of (14.171) the following holds 

HAH < C5, 1 61 77)| < C5(\x\ 2 + |f|)“. 

Therefore, from (14.211) and (14.221) it follows that 

u <b, DP >= P u>b (x,t) + w u< b(x,t), 
where P Ujb (x,t) is a polynomial of degree at most k — 1 such that 

IIA,fell < C6, \w U}b (x,t)\ < C5{\x \ 2 + | 


We also observe that in (14.201) we have that 7 0 only when \m\ + 21 < |g| + 2k < k — 1, 
and that, furthermore, the coefficient does depend on A , u and b. Now for a general polynomial 
of the form J 2 o<\m\+ 2 i<k a m/x m t e one has 

L(uP) = R(x, t ) + E iW m ,i(x,t), 


(4.23) 
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where 

(4.24) R(x,t) = ^2 d q , K x q t K , 

0<\q\+2K<k-l 

and the coefficients a m) e and d qtK of P, and satisfy (13.12p as in the case of heat equation. 
W e next introduce the appropriate notion of approximating polynomial in the present context. 

Definition 4.6. We say that a parabolic polynomial P of degree < k is approximating of order 
k at (0,0) for ^ with respect to L and g if the coefficients dn^ m of R(x,t) in the representation 
(I4.23P , (14.241) above coincide with the coefficients of the Taylor polynomial of order k — 1 for g 
at (0,0). 

With this notion in place, we now state the analogue of Lemma 13.61 

Lemma 4.7. Let u, v be as in Theorem \4.5\ Assume that for some r < 1 and P an approximating 
polynomial of order k for ^ at (0,0) with respect to L and g, with ||P|| < 1, one has 

(4-25) |\v - uP\\ L °°(Gn Qr (o,o)) < r k+1+a . 

Then, there exists an approximating polynomial P of order k such that for some C, p > 0 
universal, we have 

(4-26) \\P ~ P\\ L °°(GnQ r m) <Cr k+a , 

and 

( 4 -27) ||u - «^||Lo=(GnQ pr ( 0 , 0 )) < ( pr) k+1+a . 

Proof. The proof of Lemma 14.71 follows by arguing as in that of Lemma 13.61 We define 

v = uP + r fc+1+Q T(-, —), 

where P satisfies the hypothesis of the lemma. Since g £ H k ~ 1+a {G D ^ 2 ( 0 ,0)), from the bounds 
in (14.201) we have for all ( x,t ) £ G D Q r ( 0, 0 ), 

(4.28) \g{x,t) - P g (x,t )| < C5(\x\ 2 + | t\) L ~^ L , 

where P g is a polynomial of degree at most k — 1. 

We note that since P is approximating for — at (0,0) with respect to L and g, we have that 
P g (x,t) = R(x,t). Therefore, from (14.191) and the bounds in (14.201) we obtain in G D Q r { 0, 0) 

(4.29) | r k ~ l+a Lv(^, 4)| = | L{uP) - g\ < C5r k ~ 1+a , 

where L is as in (14.101) . The estimate (14.291) implies 

(4.30) \Lv\ < CS, 

where L is as in (14.101) . As a consequence, for a subsequence 6 —> 0 we have v = v(5) —> vo, where 
vo is as in (13.191) . Now, similarly to the proof of Lemma [3761 there exists Q 0 such that (13.201) - 
(13.231) holds. Moreover, as in the case of heat equation, the polynomial P(x, t ) + r k+a Qo(f , ^) 
need not be approximating for ^ with respect to L and g. Therefore as before, we modify Q 0 
to Q such that P(x,t) + r k+a Q(ff , is an approximating polynomial of order k for ^ at (0, 0) 
with respect to L and g. Since P is already an approximating polynomial for ^ the coefficients 
of Q should satisfy (I3.25P similarly to the situation of the heat equation. The only difference 
in the present case being that in the analogue of (13.251) the coefficients cf£ would additionally 
depend on A and b, besides u. The rest of the arguments remain the same as in the proof of 
Lemma 13.61 and the desired conclusion follows. 

□ 
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Proof of Theorem f.5. As previously, it follows by applying Lemma 14.71 repeatedly. In this 
case, in order to start the process of iteration we determine an approximating polynomial P = 
from (13.121) where d qA ’s are determined by the Taylor polynomial P g (x,t ) of 
order k — 1 for g. In view of the procedure described after (13.121) . such a polynomial can be 
determined. Then, by multiplying v and P by a suitable constant, the hypothesis of Lemma 
rm holds for some small enough universal r$. Therefore, by applying the Lemma 14.71 iteratively 
with r = ro, pro, p 2 ro and so on, we obtain a limiting polynomial Po which has the following 
representation 


(4.31) 


P 0 (x,t) = P0M) + X> 1 r 0 ) k+a Qi{ 


t 


2=1 


’ {p^rof 


where P is the above polynomial which is determined before the first step of the iteration and 
Qi are the polynomials determined after the Pth application of Lemma [4.71 Moreover, with such 
a Po, we have that (13.321) holds. The rest of the proof remains the same as that for the heat 
equation. 

□ 


5. Application to the parabolic obstacle problem 

As mentioned in the introduction, we close the paper with an application of Theorem 13.11 
to the parabolic obstacle. For a measurable set E C M n+1 , we indicate with x.q its indicator 
function. We consider the following problem studied in [ CPS ]. 

Problem 5.1. Given a domain D C M n x M, consider a function u(x,t ) defined in D such that 
u , Du are continuous, and define the coincidence set as 

A = {(x, t) € D | u(x, t ) = Du(x, t ) = 0}. 

With Ll = D\ A, suppose that u solves the equation 

A u-u t = xn- 

The free boundary is defined as 

r = r(u) = dnnD. 

With this setup, we mention a corollary of Theorem 13.11 

Corollary 5.2. Let (xo,to) € T be a point as in Theorem 13.1 and Lemma 13.3 in [CPS] , Then, 
<9C n Q\/i{xo, to) is C°°. 

Proof. By the Cl’ 1 regularity of u as in Section 4 in (CPS] , we have that the spatial derivatives 
DiU vanish on dPl for all i = 1 ,..n. Moreover, Lemma 13.3 in (CPSj implies that D^u also 
vanishes continuously at the free boundary T n Qi/a(xo, to) near (xoCo)- Now, Theorem 14.1 in 
[CPS] implies that Q im(xo, to)CT is C 1,a regular which follows by an application of the boundary 
Harnack inequality as in [ACS] . Moreover, in the proof of Theorem 14.1 in [CPS] it is evident 
that, without loss of generality, one can assume that T D Q p / 4 {xo,to) = {(x,t) \ x n = f(x',t)} 
and that the following holds, 

3 

(5.1) u n (x o + — e n , t 0 - (p/ 16) 2 ) > c 0 , 

lu 

for some Co, p > 0 universal. Moreover, p can be chosen in such a way that the point (xq + 
joe n ,to — (p/16) 2 ) is at a parabolic distance from T bounded from below by a universal constant 
Co- Now 

(5.2) 


u(x',f(x',t),t) = 0. 
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Therefore, by differentiating the equation (15.21) with respect to the variables x\,...x n -i,t, we 
obtain that 


(5.3) 


DjU 

D n u 


Dif, 


Dtu 

D n u 


D t f. 


Since (15.11) is a scaled version of the normalization (12.31) . this implies that if we take v = D^u and 
u = D n u in Theorem 12.21 we obtain from (15.31) that D'f £ H 1+a . Similarly, with v = Dtu and 
u = D n u , by application of Theorem 12.21 we find that Dtf £ H 1+a . This implies that / £ H 2+a , 
i.e., the free boundary is H 2+a regular. We now proceed inductively as follows. Suppose we 
know that / and hence the free boundary is in H k+a for some k > 2. Then, by applying Theorem 
EH to v = Dtu and u = D n u, we obtain from (15.31) that D'f £ H k+a . Similarly, with v = DfU 
and u = D n u , we find that Dtf £ H k+a . This clearly implies that / £ H k+l+a and hence the 
free boundary T D Qi/ 4 (xo,to) is H k+l+a . Therefore, we can repeatedly apply Theorem 13.11 to 
conclude that T D Qi/i{xo,tf) is smooth. 

□ 


Remark 5.3. As mentioned before in the introduction, one can in fact establish space-like real 
analyticity of the free boundary by employing the hodograph transform in [CPS] (see Theorem 
15.1 in [cps];. Nevertheless, similarly to the elliptic case, the proof of Corollary 15.21 provides 
a new perspective in the study of parabolic free boundary problems. It remains an interesting 
question to see if one can establish an analogue of Theorem 1,9.11 and Corollary 1 5. .91 in the thin 
parabolic obstacle problem studied in |PCPT| where the use of the hodograph transformation 
does not appear feasible. 
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